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The study of transport coefficients of strongly interacting matter got impetus after the discovery 
of perfect fluid ever created at ultrarelativistic heavy ion collision experiments. In this article, 
we have calculated one such coefficient viz. electrical conductivity of the quark gluon plasma 
(QGP) phase which exhibits a momentum anisotropy. Relativistic Boltzmann’s kinetic equation 
has been solved in the relaxation-time approximation to obtain the electrical conductivity. We 
have used the quasiparticle description to define the basic properties of QGP. We have compared 
our model results with the corresponding results obtained in different lattice as well as other model 
calculations. Furthermore, we extend our model to calculate the electrical conductivity at finite 
chemical potential. 
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I. INTRODUCTION 

Transport coefficients are of particular interest to quan¬ 
tify the properties of strongly interacting matter created 
at relativistic heavy ion collisions (HIC) and these 
coefficients can be instrumental to study the critical 
properties of QCD medium. The fluctuations or external 
fields cause the system to depart from its equilibrium 
and a non-equilibrium system has been created for a 
brief time. The response of the system to such type of 
fluctuations or external fields is essentially described by 
transport coefficients eg. the shear and bulk viscosities, 
the speed of sound etc. In recent years a somewhat 
surprising result in the quark gluon plasma (QGP) story 
has occurred when the practitioners in this field tried 
to satisfy the collective flow data as obtained in collider 
experiments. In order to get the required collective 
flow in the framework of viscous hydrodynamics, the 
value of shear viscosity to entropy density ratio (ji/s) 
comes out to be very small [H-Q- The tiny value of 77/s 
indicates the discovery of most perfect fluid ever created 
in laboratory. This perfect fluid is described as strongly 
interacting quark gluon plasma Thus the study 

of various transport coefficients is a powerful tool to 
really understand the behaviour of the matter produced 
in the ultra relativistic heavy ion collision (uRHIC) 
experiments at RHIC and LHC. 

Recently electrical conductivity has gained a lot of 
interest due to the strong electric field created in the 
collision zone of uRHIC experiments IM3. It has been 
observed that strong electric and magnetic fields are 
created in peripheral heavy ion collisions whose strength 
are roughly estimated as eE = eB — (where ttItt is 
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the mass of the pion) within proper time 1 — 2 fm/c [HI- 
This large electrical field can significantly affect the 
behaviour of the medium created in these collisions 
and the effect depends on the magnitude of electrical 
conductivity (cTei) of the medium. Uei is responsible 
for the production of electric current generated by 
the quarks in the early stage of the collision. The 
value of (Tel would be of fundamental importance for 
the strength of Chiral Magnetic Effect [iJI which is a 
signature of CP-violation in strong interaction. Further 
the electrical field in mass asymmetric collisions (e.g. 
Cu-Au collisions etc.) has overall a preferred direction 
and thus generating a charge asymmetric flow whose 
strength is directly related to cJei [1^. Furthermore, Ogi 
is related with the emission rate of soft photons [ij| 
accounting for their raising spectra [il,|i3. Despite of 
the importance of electrical conductivity, it has been 
studied rarely in the literature for the QGP phase. 

With the discovery of “most perfect fluid ever gen¬ 
erated” , another important observation has been made 
that this fluid possesses momentum-space anisotropies in 
the local rest frame (LRF) [13 [il. This has important 
implications for both dynamics and signatures of the 
QGP. Earlier it has been assumed a priori in the ideal 
hydrodynamics that the QGP is completely isotropic. 
However, recently dissipative hydrodynamics helps us 
to understand that the QGP created in ultrarelativistic 
heavy ion collisions has different longitudinal and 
transverse pressure. It has been shown heuristically in 
first-order Navier Stokes viscous hydrodynamics that 
the ratio of longitudinal pressure over the transverse 
pressure is: Pl/Pt = (3rT — 1677)/(3Tr -|- 8^), where 
77 = 77/s [13. Using the RHIC-like initial condition the 
value of Pl/Pt comes out equal to 0.5 and for LHG-like 
initial condition this ratio takes the value as 0.35 Hi¬ 
lt has also been shown that there exists an anisotropy 
in Pl versus Pt in second-order Israel-Stewart vis¬ 
cous hydrodynamics. Several other groups who study 
the early-time dynamics of QGD within AdS/GFT 
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framework have also shown the early-time pressure 
anisotropies and quote Pl/Pt = 0.31 or smaller BUS- 
In colour-glass condensate framework, the practitioners 
have found that the timescale for isotropization in 
classical Yang-Mills simulations is very long [2l|, [ 2 ^. 
From these observations and findings one can certainly 
assumes that the momentum anisotropy produced in 
the medium created in heavy ion collisions lasts for at 
least T <2 fm/c. Thus it is crucial to incorporate these 
momentum-space anisotropies in any phenomenological 
studies specifically for transport coefficients. 

In this article our main motivation is to calculate 
the electrical conductivity of an anisotropic QGP phase 
using the Relativistic Boltzmann’s kinetic equation. We 
have used the quasiparticle description to define the 
basic properties of QGP since earlier we have shown 
that quasiparticle description provides the proper and 
realistic thermodynamical and transport behaviour of 
QGP phase [ 2 ^ [i^. We want to provide the correct 
temperature dependence of (Xei since it is not yet 
established. Lattice calculations have obtained cTei at 
few values of temperature and these estimates vary 
widely (2^ - 1^ . Thus it is necessary to provide a proper 
temperature dependence to electrical conductivity. We 
will also extend our calculation at finite quark chemical 
potential (/iq). This is important since there is no 
guidance at finite fiq due to severe limitation of lattice 
calculations in this region. Very recently Berrehrah 
and collaborators have shown the variation of aeijT 
with respect to temperature at small but finite chemical 
potential in dynamical quasiparticle model (DQPM) [ 2 ^. 

Rest of the article is organized as follows : firstly 
we present the calculation of electrical conductivity in 
the isotropic case. In second subsection we introduce a 
momentum-space anisotropy in the distribution function 
of quarks and antiquarks and then derive the expres¬ 
sion of (Tel for this anisotropic hot and/or dense QGD 
medium using Relativistic Boltzmann’s approach. In 
the third subsection we will provide a brief introduction 
about quasiparticle model. We further explain the impor¬ 
tance of quasiparticle description of QGP in comparison 
to ideal description . Later we demonstrate results ob¬ 
tained in our model at zero and their comparison with 
the corresponding results obtained in lattice as well as in 
phenomenological calculations. We have also extended 
our model calculation at finite /i^. Finally we will give 
the summary and the conclusions drawn from this work. 


II. DESCRIPTION OF MODEL 


the longitudinal component of the electric field E : 

J = (JelE, (1) 

where the proportionality coefficient is known as electri¬ 
cal conductivity ((Tei). One can write down the expression 
for four current (J^) in a covariant form as: 

•/'" = / 9 f{x,p) -qg /(a:,p)}, (2) 

where f{x,p){f{x,p)) and q {q) are the distribution func¬ 
tion and electric charge for quark (anti-quark), respec¬ 
tively. Further g is the degeneracy factor. Let us first as¬ 
sume the case of vanishing chemical potential i.e., pq = 0. 
In this case quark and anti-quark distribution functions 
will become identical, so Eq. ([5]) takes the following form: 

J^ = ‘^<i,9f J (3) 

In the presence of infinitesimal external disturbance, the 
change in four current is : 

AJi^ = 2q^g^ J .^^^pt^6fix,p), ( 4 ) 

where Sf is the change in distribution function due to ex¬ 
ternal disturbance. One can obtain the 5f(x,p) by using 
the Relativistic Boltzmann Transport (RBT) equation, 
which is given by Hili: 

P^d^f{x,p) + qF°‘^Pfi-^f{x,p) =C[f], (5) 

where is the electromagnetic field strength tensor 
and C[f] is the collision integral, which in the relaxation¬ 
time approximation is given by : 

C[f]^-?^{f-f),{f-f)=6f (6) 

r 

where r is the relaxation time and /° is the equilibrium 
distribution function. In this approximation, RBT equa¬ 
tion becomes: 

f) 71^11 

P^df,f{x,p) + qF°‘^pi3 — f{x,p) = - ^6f ( 7 ) 

where is the fluid four velocity and in the local rest 
frame i.e., = (1,0). The equilibrium quark distribu¬ 

tion function at = 0 has the following form : 

f{x,p;T) = -^j^^, ( 8 ) 

where E = ^/p^ + rrP. 


A. Electrical Conductivity for Isotropic System 

The electric conductivity represents the response of the 
system to applied electric field. According to Ohm’s law, 
the spatial electric current (J) is directly proportional to 


Since we are interested only in electric field com¬ 
ponents of we take only two terms: = 

—E and = E. Thus the RBT equation ([7]) becomes: 


df 


df 


- q PoE • -E • p^— =- Sf, 


dp 


dp^ 


(9) 
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Q xO 

where can be solved by using the chain-rule and after 
differentiation we get the value of df as: 

= ( 10 ) 

By substituting 6f in Eq. (g]) we obtain the expression 
for (Tei: 


a 


iso _ 

el 



dPp 


X /°(1 - /°), 


( 11 ) 


where the subscript / implies summation over the 
flavors. Here we have taken up, down and strange flavors 
only. 


It is important to provide a connection between the cal¬ 
culation of electrical conductivity by thermal field theory 
and in Boltzmann’s kinetic approach. In an equilibrated 
system having volume V and temperature the zero 
frequency Green-Kubo [H, |3^ formula for the electrical 
conductivity is given by the current-current autocorrela¬ 
tion 


At finite quark-chemical potential (pg ^ 0), the distri¬ 
bution function for quark and anti-quark will be different. 
For quarks : 

= (16) 

and for antiquarks : 

f f{x,P;T, Pq) = (17) 

One can then generalize the expression for electrical con¬ 
ductivity in an isotropic medium for pq 0 as follows 


a^iPq^O) = 


X 


1 

T 


f 


d?p p^ 
(27r)^ 


k,/°(l 


(1-/T) -(18) 


B. Electrical Conductivity for Anisotropic System 


^el — 


pOO 

13V / {ji{0)j^{t))dt, 

Jo 


( 12 ) 


where the repeated spatial-indices, i do not imply sum¬ 
mation. In the local rest frame of fluid, the electric cur¬ 
rent density can be read as : 


j\t) = 



where M is the number of particle-species and Nk is the 
number of particles of A:-th species. For all equivalent di¬ 
rections (/), ji {t) was obtained from Boltzmann Approach 
to Multi-parton Scatterings and then the autocorrelation 
function, C{t) of the electric current density in equilib¬ 
rium has been extracted 0. For example, the variance, 
(7(0) has been computed analytically 


1 ^ 

C{Q) = ^^qlnk, (14) 

^ k=l 

which finally gives the expression for (Tgi • 

1 “ 

^ X! , (15) 

k=l 


The introduction of an anisotropic distribution func¬ 
tion is needed to properly describe the QGP created 
in heavy ion collisions [^, |^. Partons are produced 
from the incoming colliding nuclei just after the collision 
at proper time r = tq « it can be assumed that 

the newly produced partons follow an isotropic (but 
not necessary an equilibrium) momentum distribution. 
Here Qs is the gluon saturation scale. The early-time 
physics is mainly governed by the hard gluons with 
the momentum at the saturation scale which have very 
large occupation numbers of order 1/as (c^s << 1) |36l - 
[ 3 ^. For r > Q~^, the hard gluons would follow the 
straight-line trajectories and isolate themselves in beam 
direction as if no any interaction exists. As a result, the 
longitudinal expansion causes the medium to become 
much colder in the longitudinal direction than the 
transverse direction, ie. p± » Pz ^ 1 /t and a local 
momentum anisotropy appears. 

The anisotropic distribution can be obtained by 
stretching or squeezing an isotropic one along a certain 
direction, thereby preserving a cylindrical symmetry in 
momentum space. In particular, the anisotropic distribu¬ 
tion relevant for HICs can be approximated by removing 
particles with the large momentum component along the 
direction of anisotropy, n as [1^ : 


which is nothing but the non-relativistic Drude’s for¬ 
mula to calculate the (Je\- Similarly one can show that 
the expression for del obtained by solving relativistic 
Boltzmann kinetic equation (see Eq. (11)) can also 
be approximated as Drude’s formula assuming small 
electric field and no cross effects between heat and 
electrical conductivity [ 13 . 


/aniso(p) = fiso + C(P ' n)^^ , (19) 


where fiso is an arbitrary isotropic distribution function 
and f is the anisotropic parameter and is generically de¬ 
fined as: 


(p^t) 

2(pi) 


- 1 , 


( 20 ) 
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where = P-n and Pt = P — n(p.n) are the com¬ 
ponents of momentum parallel and perpendicular to n, 
respectively. There have been significant advances in the 
dynamical models used to simulate plasma evolution with 
the momentum-space anisotropies [4ll - l4^ . Recently two 
of us have investigated the effects of anisotropy on the 
quarkonia states by the leading-anisotropic correction to 
the potential at T=0 jH, 

If fiso is a thermal ideal-gas distribution and ^ is small 
then ^ is also related to the shear viscosity of the medium 
via one-dimensional Bjorken expansion in the Navier- 
Stokes limit 


medium, after summing over the flavours (/) 



= (26) 


where the anisotropic term, A is given by 


10 r] 
Tt s ’ 


( 21 ) 


In an expanding system, non-vanishing viscosity im¬ 
plies finite-momentum relaxation rate and therefore 
an anisotropy of the particle momenta appears. For 
rjjs = 0.1 — 0.3 and tT =1 — 3, one finds that ^ = 1. 


/ 

Like in the isotropic medium, we can generalize the 
electrical conductivity in an anisotropic medium for pq 
0 ; 


As we have explained, hot QCD medium due to expan¬ 
sion and non zero viscosity, exhibits a local anisotropy in 
momentum space, and the quark distribution function 
(or Fermi-Dirac distribution function) takes the follow¬ 
ing form for fiq = 0: 


/aniso(()05 Pi F) — 


j(-\/p^+C(P-n)^+™^)/r _|_ I 


( 22 ) 


For weakly anisotropic systems << 1), one can expand 
the quark distribution function as follows : 


/aniso((^i Pi F) — 


1 


+ 1 

0 ? 


2ET 


(p.n) 


oE/T 


(g£/T + l)2’ 


= r- 


2ET 


(P • n) 


(23) 


where p = (psin0cos()),psin0sin(/),pcos0) and n = 
(sin a, 0, cos a), a is the angle between p and n. Af¬ 
ter substituting the anisotropic distribution function in 
Eq. (fTTl) . the expression of electrical conductivity in 
anisotropic medium system is modified as: 

= 5 /^ / U^soil - faniso) ( 24 ) 

On neglecting the higher-order distribution function 
(e.g., /o etc.) we get the expression for electrical con¬ 
ductivity as: 


^aniso ^ 

C^el = 9f 


- ^9f 


T 

9/ 


(Pp p^ 
d?p p^ 

r2 J (27r)3 


r//°(l-/°) 


(25) 


Jif) 


0^2 E/T 


{p-nf 


Now using the definition of p and n, and integrating over 
0 and 4>, we get the electrical conductivity in anisotropic 


atr°{pq^0) = 


2tt‘^T 


! dp^ 


f 


r,/“(l-/“)+r,/“(l-/7) 
,6 


^Y^Y.9f9^f [dp^ 


(28) 


/ 






where and acquire the form as given in Eq. (16) 
and ( 17 ), respectively. 


C. Quasiparticle Model: Effective Masses and 
Relaxation Times 


Quasiparticles are the thermal excitations of the inter¬ 
acting quarks and gluons retaining the quantum numbers 
of the real particles, i.e., the quarks and gluons. In quasi¬ 
particle models [i^, QGP is described by the system of 
massive noninteracting quasiparticles where the mass of 
these quasiparticles is temperature-dependent and arises 
because of the interactions of quarks and gluons with the 
surrounding matter in the medium. The effective mass 
of these quasiparticles is given by [23|: 


A ^ (29) 

where mo,/ is the current mass of the flavour, / and mth.f 
is the thermal mass of the flavour, /, which is given by 


^thj - - 7 - 



Al] 

7r2T2 J ■ 


(30) 
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Here is the QCD running coupling constant which in 
two-loop has following form [SOl l5ll | : 


OisiT, flq) = 


X 


g'^jT, gq) ^_for_ 

(33 - 2Nf) In 


where Kt is the QCD scale-fixing parameter which char¬ 
acterizes the strength of the interaction. It originates 
from the lowest non-zero Matsubara modes [l^. Here 
parameter a is equal to 


In Eqs. (fTSl) and (I^Hl) . t/ is the collision time for which 
we use the following expressions for quarks (anti-quarks) 
from Ref. 1^: 


1 


^ 9 ( 9 ) 


IbalTlog (^) (l + 0.06iV( 


eff 

F 


(32) 


where A^|? is the number of effective flavour degrees of 



FIG. 1: Variation of QPM and ideal distribution function with 
respect to momentum at two different temperatures. Solid 
and dash-dotted curve represents the light quark distribution 
function in QPM at T — 180 and 500 MeV, respectively. 
Similarly short-dashed and long-dashed curve shows the cor¬ 
responding results in ideal description. 

freedom. 

Before going to the electrical conductivity results, we 



FIG. 2: Variation of QPM and ideal distribution functions 
with respect to temperature at two different momentum.Solid 
and dash-dotted curve represents the light quark distribution 
function in QPM at p = 0.01 and 1000 MeV, respectively. 
Similarly short-dashed and long-dashed curve shows the cor¬ 
responding results in ideal description. 


want to just provide the difference between the ideal de¬ 
scription (where the mass consists of only the current 
mass) and the quasiparticle description (where the mass 
consists of current mass along-with a thermally gener¬ 
ated mass) of QGP. To understand the crucial difference 
between these two descriptions, it is better to plot the 
occupation probability or distribution function in both 
descriptions. In Fig. I, we have shown the variation 
of QPM distribution function for light quarks {u and/or 
d) with respect to momentum at two different tempera¬ 
tures T = 180 and 500 MeV. We further compare these 
results with the corresponding ideal distribution func¬ 
tions where we use only the current mass of the quarks 
in the mass term and there is no any thermal contribu¬ 
tion. We observe that at low momentum the difference 
between QPM and ideal description is large. The QPM 
occupation probability is low in comparison to ideal case 
and thus we can indirectly say that the number density 
will also be small. However, at higher momentum, both 
picture give the same distribution function. Similarly, 
Fig. 2 demonstrates the variation of QPM distribution 
function with respect to temperature at two different mo¬ 
mentum p = 0.01 MeV (low momentum) and 1000 MeV 
(high momentum). We have also plotted the correspond¬ 
ing ideal distribution function for comparison. Here we 
observed that the QPM distribution function or QPM 
occupation probability is small in comparison to ideal 
case at low momentum and this difference increases as 
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we move from higher temperatures towards lower tem¬ 
perature. However, at higher momentum the difference 
between QPM and ideal distribution function is quite in¬ 
significant over the entire temperatures range considered 
by us. As we know from non-relativistic Drude’s formula 
that the electrical conductivity is directly proportional 
to the number density which is nothing but the integra¬ 
tion of distribution function over momentum space at 
any particular temperature. Therefore one can under¬ 
stand the usefulness of QPM description in describing 
the thermodynamical and transport properties of QGP 
specially near the critical temperature (Tk). 

III. RESULTS AND DISCUSSIONS 



FIG. 3: Variation of ratio of electrical conductivity to tem¬ 
perature ((j|f‘®°/r) with respect to temperature at different 
values of anisotropy parameter ^ = 0.0, 0.3 and 0.6 and they 
are shown in the plot by solid, dashed and dash-dotted line, 
respectively. Green points are the results as obtained from dy¬ 
namical quasiparticle model (DQPM) in Ref. [^. All other 
symbols are the data points obtained from various lattice cal¬ 
culations [2 r1 - [^ . 

In Fig. 3, we present the variation of ratio of electrical 
conductivity to temperature (o'g“®°/T) with respect to 
temperature at p,, = 0 for an anisotropic QGP. For 
weakly anisotropic system, we choose three constant 
values of ^ as 0.0, 0.3 and 0.6, where ^ = 0.0 represents 
tJei/T in isotropic case as mentioned in Eq. (If). We 
have compared our model results with the corresponding 
results obtained in various lattice calculations [2^ -1 ^ . 
We further compared our model results with the aei/T 
obtained in DQPM model Q- We observe that /T 


increases monotonically with temperature starting from 
a lower value at T = Tc, where Tc is the crossover 
temperature for the transition from QGP to hadron 
gas (HG). Further we observe that as the anisotropy 
increases from 0.0 to 0.6, the value of a^xjT decreases for 
all the temperatures. (Tg“®°/T satisfies the lattice results 
well when the anisotropy parameter has a value equal 
to 0.6. The results obtained from DQPM overestimate 
the value of u^xjT as compared to the lattice as well as 
our model calculations. However, we cannot say at this 
moment the exact status of any model calculations since 
the lattice results are distributed over a wide range. 



FIG. 4: Variation of ratio of electrical conductivity to temper¬ 
ature with respect to temperature at two different 

values of anisotropy parameter ^ = 0.0 and 0.6 for quark 
chemical potential = 0 and = 200 MeV. 


Fig.4 represents the variation of with 

respect to temperature at quark chemical potential 
jjLq = 200 MeV and compare them with our model 
results at fj,q = 0 MeV. Here we have shown the 
results for two different values of ^ equal to 0.0 and 0.6. 
Solid and long-dashed curve demonstrate the at 

fiq = 0.0 for ^ = 0.0 and 0.6, respectively. Short-dashed 
and dash-dotted curve show the results at ^q = 200 MeV 
for ^ = 0.0 and 0.6, respectively. At large temperatures 
we found that at finite ^q remains same as 

in the case of zero /r,. However, the value of 
becomes large for finite /ig in comparison to the value 
at = 0 MeV as the temperature decreases below 
T = ATc- This behaviour is quite consistent with the 
results obtained in a very recent DQPM calculations at 
finite /ig [^ . 
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In summary we have studied the behaviour of electri¬ 
cal conductivity for QGP phase in the framework of rel¬ 
ativistic Boltzmann’s kinetic equation using relaxation 
time approximation. In this work we describe the QGP 
as a system of quasiparticles having temperature and 
chemical potential dependent mass along with their rest 
masses. Firstly we derive the expression to calculate del 
in an isotropic medium. Later a momentum anisotropy 
in the distribution function of quarks and anti-quarks has 
been introduced and taking the leading order contribu¬ 
tion in distribution function we have obtained del in an 
anisotropic medium which is equal to the dg[° minus the 
correction factor due to momentum anisotropy. We have 
plotted dg“®°/T for different values of ^ which is equal 
to 0.0, 0.3 and 0.6. We have compared our model results 
with the corresponding results obtained in various lat¬ 
tice results and found a reasonable agreement between 
them for certain value of Further, we have plotted 
our model expectations for d“‘®°/r at /r, = 200 MeV 
and compare them with the model results obtained at 
zero chemical potential. Both in isotropic as well as in 
anisotropic case, we find the similar behaviour as ob¬ 


served in DQPM calculations i.e., increase in the 
value of dei/T with increase in quark chemical potential 
near Tc- At very high temperatures (T ~ 8 — lOTc), the 
difference between /r, = 0.0 and fig = 200 MeV case is 
very small. These calculations are done in a static system 
when there is no any proper time dependence has been 
given to the anisotropy parameter. However, in realis¬ 
tic situation, ^ varies with the proper time starting from 
the initial proper time up to a time when the system be¬ 
comes isotropic and ^ becomes zero. Thus one has to 
incorporate a proper time dependence to the anisotropy 
parameter. Work in this direction is in progress and will 
be presented elsewhere. 
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